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CURRENTS  ON  GENERALIZED  YAGI  STRUCTURES 


INTRODUCTION 

This  report  evaluates  properties  of  currents  on  an  infinite  linear  array  of  parasitic  radia- 
tors when  one  of  the  radiators  is  excited.  In  practical  applications  Yagi  antennas  use  arrays 
of  parasitic  electric  dipoles.  Therefore  the  detailed  analysis  deals  with  a generalized  Yagi 
(dipole)  structure.  The  structure  has  been  generalized  by  allowing  the  dipole  antennas  to  be 
inclined  at  an  arbitrary  angle  to  the  array  axis;  in  a conventional  Yagi  the  dipoles  are  per- 
pendicular to  the  array  axis. 

When  the  array  is  excited,  the  currents  at  each  radiator  may  be  decomposed  into  a 
surface-wave  component  plus  a correction  component.  The  surface-wave  currents  are  of 
particular  importance  in  design,  and  their  wave  properties  are  found  for  various  angles  of 
inclination  of  the  dipoles. 

Calculations  are  carried  out  within  the  framework  of  a network  formulation.  The  in- 
tegral solution  obtained  from  this  formulation  is  used  to  demonstrate  the  interrelation 
between  the  properties  of  the  array  of  short-circuited  radiators  excited  parasitically  and 
those  of  the  same  structure  when  each  radiator  is  excited  by  a real  generator,  as  in  a phased 
array.  Analytical  (closed)  forms  previously  obtained  in  phased -array  studies  are  used  to  elimi- 
nate numerical  difficulties  due  to  slow  convergence  of  the  series  which  arise  in  previous 
treatments  of  long  Yagi  antennas.  For  the  special  case  of  dipoles  inclined  at  the  angle 
arctan-v/T13  aresin  >^2/3  a 54.74°  to  the  array  axis,  the  functional  form  of  the  solution  sim- 
plifies remarkably.  The  surface-wave  and  (feed)  correction  currents  are  then  explicitly  eval- 
uated, in  the  complex  plane,  as  a pole-residue  contribution  and  a branch-cut  contour  integral. 

As  recounted  by  Professor  Uda  1 1,2) , the  Yagi-Uda  antenna  was  invented  in  1926. 
Further  practical  and  theoretical  studies  were  undertaken,  but,  as  noted  by  Ehrenspeek  and 
Poehler  (3) , in  tire  late  1950’s  there  existed  no  rigorous  solution  of  the  Yagi  problem.  The 
experimental  results  were  restricted  to  special  eases,  with  no  attempt  made  to  find  a connec- 
tion between  them.  Ehrenspeek  and  Poehler  developed  genera!  design  principles  for  long 
Yagi  antennas.  Their  experiments  demonstrated  the  dominant  role  played  by  the  surface- 
wave  parameters  in  determining  the  performance  of  this  antenna.  This  is  now  well  understood 
within  the  context  of  surface-wave  antenna  design  (4) . 

The  variation  of  the  phase  velocity  of  the  surface  wave  on  infinite  Yagi  structures  was 
analyzed  by  Scngupta  (5).  Mailloux  [6,7)  provided  a complete  solution,  including  excita- 
tion coefficient,  for  the  infinite  Yagi  structure  excited  at  one  element  of  the  structure.  He 
then  applied  these  results  to  finite  Yagi  arrays,  obtaining  excellent  agreement  with  experi- 
ment and  with  an  alternative  theory  for  such  arrays  by  King  and  Sandler  [ 8) . Gatcly  et  al. 

[9)  showed  that  a comparatively  simple  calculation  which  treats  the  dipoles  of  the  Vagi 
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array  as  minimum -scattering  antennas  [10,11]  also  yields  similar  results  for  finite  arrays.  For 
any  given  array  configuration,  even  if  it  includes  a large  number  of  elements,  such  direct 
calculations  of  radiation  from  a finite  array  no  doubt  now  provide  the  most  convenient 
route  to  accurate  results.  On  the  other  hand,  because  of  the  many  variables,  estimates  for 
appropriate  designs  continue  to  be  most  readily  derived  from  the  surface-wave  point  of  view. 

In  the  next  section  of  this  report  the  problem  posed  by  a linear  array  of  identical 
parasitic  radiators  is  given  a generic  network  formulation.  A formal  solution  in  the  form  of 
an  integral  is  obtained  for  the  currents  produced  on  such  a structure  when  a single  element 
is  excited  (the  Green’s  function  currents.)  The  technique  for  solution  was  described  by 
Mailloux  [6]  and  follows  closely  that  employed  in  the  analysis  of  phased  arrays  [12,13] . 

In  the  present  report  it  is  shown  that  the  active  impedance  (the  input  impedance  to  any 
antenna  element  when  all  are  excited  with  uniform  amplitude  and  uniform  phase  difference  £) 
is  intimately  connected  with  the  surface-wave  parameters.  In  particular,  surface  waves  can 
occur  only  for  values  of  £ denoted  £s  corresponding  to  the  “invisible  region”  in  phased-array 
parlance.  For  such  values  no  radiation  occurs,  and  the  active  impedance  is  reactive.  The  ratio 
of  surface-wave  velocity  to  the  velocity  of  light  is  kD/£s,  where  kD  is  the  element  spacing  in 
electrical  radians  and  £s  is  a zero  of  an  active  impedance  as  a function  of  phasing  angle.  It 
would  seem  that  analyses  carried  out  on  a variety  of  phased-array  structures  can  now  be 
turned  to  account  in  the  design  of  surface-wave  antennas. 

In  Refs.  12  and  13  Wasylkiwskyj  and  Kahn  analyzed  an  infinite  linear  array  consisting 
of  dipoles  (minimum-scattering  antennas)  oriented  at  an  arbitrary  angle  with  respect  to  the 
array  axis.  An  essentially  closed  form  for  the  active  impedance  is  available  from  Ref.  12. 

This  form  simplifies  remarkably  for  dipoles  inclined  at  the  angle  Oq  a arcsin\/2/3  = arctan>/2~ 
The  special  nature  of  this  angle  for  dipoles  was  first  noted  by  Hazeltine  [ 14]  .*  The  third 
section  of  this  report  is  devoted  to  the  Yagi  structure  consisting  of  dipoles  inclined  at  this 
special  angle.  For  this  case  tire  formal  integral  solution  can  he  evaluated  by  contour  integra- 
tion in  the  complex  plane.  The  surface-wave  components  of  currents  at  the  antenna  terminals 
are  evaluated  as  pole-residue  contributions,  and  the  remaining  components  of  current  are 
evaluated  as  branch-cut  integrals.  On  the  infinite  structure  these  branch-cut  “correction"  or 
space-wave  components  of  current  are  the  only  ones  giving  rise  to  radiation  away  from  the 
structure.  In  contrast  to  the  surface-wave  components,  which  retain  a constant  amplitude 
along  the  structure,  the  correction  components  decrease  with  distance  for  antenna  elements 
removed  from  the  one  excited.  It  is  found  that  this  decrease  is  at  least  as  rapid  as  1 In,  where 
n is  number  of  elements  removed  from  the  excited  element.  In  a Yagi  array  Um  type  of  cur- 
rent would  be  responsible  for  “feed  radiation." 


•Hareltine  (14 1 discovered  that  the  italic  coupling  between  parallel  magnetic  cell*  (more  generally  dipole*) 
wa*  eliminated  when  the  coila  are  inclined  to  a line  connecting  their  center*  at  the  special  angle  6q  • 
sretan  s/2.  Not  only  the  static  coupling  term  inversely  proportional  to  the  cube  of  the  coil  separation  but 
also  a dynamic  term  inversely  proportional  to  the  square  of  the  separation  are  eliminated.  Only  the  dynamic 
mdiation  term  inversely  proportional  to  the  separation  remains.  The  author  is  grateful  to  Or.  Harold  A. 
Wheeler  for  bringing  Haaeltine's  discovery  to  his  attention.  (Or.  Wheeler  has  discussed  Hardline ’•  discovery 
in  a recent  book  Haiti  line  the  Profetsor,  published  by  the  Hardline  Corporation,  Greentawn,  NY.,  in 
197ft.) 
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The  last  section  deals  with  the  parameters  of  surface-wave  components  on  the  general- 
ized dipole  Yagi  structures.  On  the  one  hand,  except  at  the  one  special  angle  treated  in  the 
third  section,  the  complexity  of  the  integrand  in  the  formal  solution  would  seem  to  preclude 
a corresponding  complete  evaluation.  On  the  other  hand  the  form  obtained  for  the  active 
impedance  [12]  is  convenient  for  computation.  Consequently  it  becomes  economical  to 
compute  the  actual  variation  of  surface- wave  parameters  with  frequency  for  any  particular 
structure  of  interest.  Examples  of  such  calculations  are  given.  For  conventional  Yagi  struc- 
tures (0g  = 90°)  surface-wave  solutions  have  been  found  only  for  capacitive  (short)  dipoles. 

It  is  shown  that,  for  sufficiently  small  angle  of  inclination  9q,  surface  waves  can  also  exist  on 
arrays  of  inductive  dipoles.  In  the  computations  the  mutual  coupling  between  antennas  is 
(apart  from  a scale  factor)  approximated  as  the  coupling  between  minimum  scattering 
antennas  having  radiation  patterns  of  short  dipoles.  Therefore  results  obtained  for  inductive 
dipoles  apply,  strictly  speaking,  to  short  inductively  loaded  dipoles.  Because  of  the  slow 
change  of  the  radiation  pattern,  the  theory  of  mutual  coupling  between  minimum  scattering 
antennas  indicates  that  these  results  approximately  apply  to  arrays  of  (unloaded)  dipoles 
somewhat  longer  than  1/2  wavelength. 


NETWORK  FORMULATION 

An  infinite  array  of  identical  antennas  uniformly  spaced  along  a straight  line  or  axis  is 
characterized  at  the  terminals  of  the  antennas  by  voltages  Vm  and  currents  /„,  where  the 
subscripts  denote  the  mth  and  nth  antennas.  Those  terminal  quantities  are  related  by  the 
open-circuit  impedance  coefficients  Zmn ; 

i - £ z„„  /„ . (i) 

n«-«* 

In  view  of  the  symmetry  of  this  array,  the  impedance  coefficients  depend  only  on  the  separa- 
tion between  the  mth  and  nth  antennas  (n  - m)V,  where  D is  tire  separation  of  adjacent 
antennas.  Thus  Zmn  is  a function  of  only  the  difference 

e » n - m.  (2) 

Lo rents  reciprocity  in  the  electromagnetic  field  implies  that  the  matrix  of  impedance  coeffi- 
cients is  symmetrical: 


It  follows  that 


Zmn  “ %nm 


Zmn  = Zv  a Z_v 


The  basic  network  relation  (1)  can  therefore  be  rewritten  as 

35  S , Zy  Im*v  ■ 

v»— 

This  is  an  infinite-order,  linear,  finite-difference  equation. 


(3) 

(4) 

(5) 
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Let  us  assume  that  only  one  of  the  antennas,  the  £th  antenna,  is  excited  by  a voltage 
generator  Eg  with  internal  reactance  Xa  (Fig.  1).  The  remaining  antennas  are  each  terminated 
in  a like  reactance  Xa ; they  are  excited  only  parasitically  through  mutual  coupling.  These 
conditions  constrain  the  terminal  voltages  and  currents  to 

~ Eg  ^m2  ~ J Xa  Im  , (6) 

where  6mg  is  the  Kronecker  delta: 

sm£  = 1.  m = £,  (7a) 

=0,  m * fi.  (7b) 

When  these  constraints  are  inserted  into  the  difference  equation  (5),  the  result  is 

Eg  ^ (Zy  + J Xa  &0i>  (8) 

To  solve  the  difference  equation  (8),  one  can  introduce  the  Fourier  transform  of  the 
currents  Iv : 

« 

■*«)  - £ (9a) 

and,  inversely, 

(9b) 


Fig.  I — Infinite  *my  of  identic*!  tntecuu*, 
with  one  inUtuu  excited 
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From  this  definition  it  follows  that  the  Fourier  transform  of  the  currents  Im  + is  e_;m 


Im+v 


ime'M]  dt 


(10) 


The  Kronecker  delta  is  given  by 


$m£  = 


e-Um-m 


dt. 


(11) 


On  insertion  of  the  integral  representations  (10)  and  (11),  difference  equation  (8)  becomes 


i f d*  = E + Wo m e~/(m+y)*  d$-  <12> 

Interchanging  the  order  of  summation  and  integration  and  equating  the  resultant  transforms, 
one  may  solve  for  => 

m Ee*m 

■ Txf^m  • (13> 

where 

Q(0  B • d4) 


The  superscript  £ makes  explicit  the  dependence  of  the  solution  on  excitation  at  the  |K»rt  of 
the  £th  antenna  only.  The  currents  /^)  which  constitute  the  solution  of  (8)  (in  effect  a 
Green's  function  solution)  may  then  be  recovered  from  (13): 


/<«) 

n 


2*]m  jXc  ♦ Q($) 


(16> 


It  is  clear  from  (15)  that  the  properties  which  distinguish  tire  terminal  currents  on  one  parti- 
cular  array  from  those  on  another  reside  entirely  in  Xa  and  the  functional  form  of  Q(|)  (15] . 
From  symmetry  considerations  or  from  (15) 


W • «l  ■ 


(16) 


In  the  test  of  this  report  the  superscript  6 will  be  suppressed  by  setting  t°0. 

In  the  phased- array  literature  (12,16,17)  <?(£)  is  identified  as  tire  active  impedance,  which 
is  the  input  impedance  at  any  element  of  the  infinite  array  when  all  elements  are  excited  by 
generators  of  uniform  amplitude  and  constant  phase  difference  (£  radians).  Therefore  an 


5 


KAHN 


i-  ..nate  connection  exists  between  the  characteristics  of  the  parasitic  array  and  those  of  the 
same  array  of  antennas  employed  as  an  active  phased  array.  The  implications  of  this  connec- 
tion will  be  developed  subsequently. 

A significant  difficulty  arises  in  the  evaluation  of  the  currents  In  from  (15)  when  the  in- 
tegrand has  singularities  for  values  of  £ in  the  interval  -7r  < £ < tt.  As  will  be  discussed  from 
several  points  of  view,  particular  singularities  contribute  surface-wave  components  of  current. 
These  particular  singularities  occur  at  the  roots  of 

}Xa  + Q(£)  = 0.  (17) 

Clearly  a necessary  condition  for  occurrence  of  a root  in  (17)  is  that  the  real  part  of  the  active 
impedance  Q(£)  vanish: 


Re  Q(£)  =>  0.  (18a) 

It  is  well  known  that  for  closely  spaced  arrays,  Q<kD<ir,  symmetry  considerations  dictate 
a range  of  values  of  $ in  which  (18a)  holds,  the  so-called  “invisible"  region  of  phasing  angles 
[12,16,17] 


kD<\(  Kir.  (18b) 

Therefore  real  roots  of  (17)  quite  generally  He  in  that  range.  Further,  it  can  be  shown  that, 
due  to  the  reciprocity  conditions  (3)  and  (4),  these  roots  occur  in  pairs  with  opposite  sign 
(Appendix).  If  ^ is  a real  root  of  (17),  then  is  also  a root,  which  will  be  denoted 

$(-P>  . p » tl.  *2 (19) 

Surface-wave  components  of  current  (components  that  maintain  a constant  magnitude 
from  one  antenna  element  to  the  next  and  a constant  phase  difference  from  one  dement  to 
the  next)  are  found  on  closely  spaced  arrays.  These  waves  have  phase  velocity 


. “D.  « CM 

* $(P)  $(P)  ’ 


(20) 


where  c is  the  velocity  of  light.  In  view  of  (18b)  they  are  “slow  waves." 

The  total  currents  /„  are  conveniently  dissected  into  surface-wave  and  correction 
components 


ln  “ • (21») 

where 

■ ',0  <*»M 

/ is  by  deftnition  a surface-wave  component  of  current  (with  the  superscript  p being 
omitted,  because  only  one  surface  wave  is  assumed  for  simplicity).  The  absolute  value  in 
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(21b)  is  justified  by  the  symmetry  of  the  problem.  The  unknown  amplitude  7sq  can  be  de- 
termined from  the  Fourier  transform  of  (21a) 


•?($)  = •*,«)  + 


where 


1 _ e« 


- + 


(22a) 


(22b) 


with  the  requirement  that  ,?c(£)  be  nonsingular  at  (•  * $ [6] . If  no  other  (type  of)  singularities 
are  present,  the  currents  / may  then  be  evaluated  straightforwardly  from  .4>#(()  and  the 
inverse  transform  (9b).  Wfien  other  types  of  singularities  do  occur,  these  must  be  examined. 

If  they  are  integrate,  no  further  analytical  difficulties  (but  quite  possibly  further  numerical 
difficulties)  may  be  involved  in  the  evaluation  of  7cr|: 

7Cfl  = ~ ) e-'"«  d{.  (22c) 


An  alternative  means  for  evaluation  of  the  currents  is  accessible  when  the  analytical 
properties  of  {?(£)  can  be  determined  in  detail.  Extension  of  <?(£)  into  the  complex  f plane 

f « W1.  > 0,  (23) 

brings  with  it  the  elegance  and  power  of  function-theoretic  techniques  for  integration.  As  a 
function  of  the  complex  variable  f,  the  active  impedance  function  will  be  denoted  by  q($). 
When  q * 1, 


q(S)  * q(nS)  - <?(*).  (24) 

The  currents  In  are  then  given  by  the  contour-integral  formula 

i Af  r("M) 

f«  " 2*7  J /*,  ♦ • (25) 

c 

which  is  equivalent  to  ( 16)  when  J is  the  unit  circle  in  the  f plane,  property  indented  about 
any  singularities.  The  proper  indentations  may  be  determined  by  introducing  a small 
amount  of  dissipation  (which  moves  the  singularities  off  the  unit  circle)  and  then  pasting  to 
the  limit  of  zero  loss  (4) . An  example  of  such  a contour  will  be  shown  in  Fig.  4a  in  connec- 
tion with  the  particular  case  of  inclined  dipoles  considered  in  the  next  sec  linn. 
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When  the  integral  (25)  is  evaluated  by  deforming  the  contour  C about  the  singularities 
of  the  integrand,  the  residue  contribution,  due  to  a simple  pole  at  f on  the  unit  circle, 
is  readily  identified  as  a surface-wave  current.  In  general 


* sn 


t-t, 


- }E, 


(26a) 


(26b) 


The  preceding  formulation  is  readily  generalized  to  an  infinite  array  of  a given  complex 
of  antennas  [13) . The  given  complex,  termed  a subarray  of  antennas,  may  comprise  one  or 
more  multiport  antennas  (Fig.  2)  or  simply  an  ordinary  subarray  of  several  identical  one-port 
antennas.  Suppose  the  antennas  in  each  subarray  have  M input  ports.  Then  the  impedance 
coefficients  of  the  infinite  structure  may  be  denoted  by  2™  n , 7n,  7i  »>  1 , 2, . . . , Af,  when?, 
for  example,  m identifies  the  with  port  of  the  mth  subarray.  The  Lorente  reciprocity  condi- 
tion (3)  is  then  replaced  by 


'in'Sunn 


^ n Mi «« 


(27a) 


Fig.  2 — Array  of  *ubamy 
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With  the  adoption  of  matrix  notation  within  the  subarray, 


Z = 1 , (27b) 

~mn  —nm 


where  Z denotes  the  transpose  of  the  M-by-M  matrix  Z . From  symmetry  (2)  it  follows 
that  ~nm  ~nm 


z = z = z . 

—mn  —v  —-v 
Accordingly  the  basic  network  relation  can  be  written  as 

where  V and  I are  Af -dimensional  column  matrices: 

— m —n 


V = 


= >U 


m 2 ;m 


Af.TJ 


and 


(28) 


(29) 


(30a) 


'2;n> 


rAf;*i  ' 


(30b) 


The  solution  of  (29)  proceeds  in  a fashion  entirely  parallel  to  the  solution  of  (5).  In  particular 


(31) 


is  obtained  parallel  with  (15).  Surface-wave  currents  may  occur  and  are  again  attributable  to 
singularities  of  the  integrand  in  (31),  that  is,  to  real  roots  of  the  determinant: 


det 


+ £(*)] 


a 0. 


(32) 


APPLICATION  TO  AN  INFINITE  ARRAY  OF  INCLINED  DIPOLES 

An  infinite,  uniform,  linear  array  of  electric  dipoles,  each  inclined  at  angle  0Q  to  the 
array  axis,  is  shown  in  Fig.  3a.  Let  the  input  impedance  to  one  of  the  dipole  antennas  taken 
from  the  array  be 
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Fig.  3a  — Infinite  linear  array  of  metal  rods  (radius  a and 
length  b)  or  short-circuited  electric  dipole  antennas 


equivalent  actual 

ANTENNA  TERMINALS 

CMS  ANTENNA 

TERMINALS  ACTUAL  ANTENNA 

I TERMINALS 


ACTUAL  LOAD 
(SHORT  CIRCUIT) 


v - 

EQUIVALENT  LOAD 
FOR  CMS  ANTENNA 


Fig.  3b  — Dipole  and  equivalent  canonical-minimum- 
scattering  (CMS)  antenna  circuit 


K * I*. 


(33) 


in  free  space.  To  use  theory  developed  for  minimum-scattering  antennas  to  calculate  mutual 
impedance  (9,11,17) , each  antenna  is  modified  through  the  addition  of  the  series  reactance 
circuit  shown  in  Fig.  3b.  At  terminal  bb'  each  dipole  Is  closely  modeled  by  a canonical- 
minimum -scattering  (CMS)  antenna  with  input  impedance  Rfl.  Terminals  cc'  are  in  all  respects 
equivalent  to  the  original  antenna  terminals  oa‘.  With  respect  to  the  terminals  bb’  the  im- 
pedance coefficients  for  the  array,  normalized  to  Rfl,  are  given  by  (11,12) 


h(2) 


sin20ft  + 3~h 


h(2)  (feDiPi) 


cos20o,  *>^0, 


“ 1,  M“0, 

whore  the  particular  spherical  Hankel  hmetions  are  (18) 


(34a) 

(34b) 
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and 


(2)  = j 


(35a) 


(35b) 


and  where  k = 2 tt/X,  D is  the  distance  between  adjacent  dipole  centers,  and  D I v I = D I n - m I 
is  the  distance  between  the  mth  and  nth  dipole  centers.  It  was  noted  following  (1)  that  the 
mutual  impedance  depends  only  on  the  distance  between  the  mth  and  nth  dipole  centers  and 
thus  on  kD  I m - n I = kD  I v I.  Impedance  quantities  without  the  caret  are  normalised  to  Rfl. 

An  infinite  Yagi  structure  is  obtained  when  the  dipoles  are  short-circuited  at  the  termi- 
nals aa  of  the  antennas.  Putting  the  short  circuit  at  the  terminals  cc  shows  that  the  Yagi  struc- 
ture is  equivalent  to  the  array  of  CMS  antennas  at  terminals  bb‘  terminated  in  reactances 
When  an  ideal  generator  is  placed  in  series  with  this  reactance  at  one  of  the  antennas,  the 
problem  of  finding  the  terminal  currents  at  all  antennas  is  a particular  example  of  the  network 
problem  formally  solved  in  the  preceding  section.  As  was  pointed  out  in  that  section,  the 
features  distinguishing  various  infinite  Yagi  structures  are  concentrated  in  the  functional  form 
of  the  active  impedance  Q(£)  for  all  phasing  angles  £,  -rr  < £ < it,  giver,  by  Eq  , (14).  A highly 
convergent  form  of  this  impedance  was  previously  obtained  by  Wasylkiwskyj  and  Kahn  [12] . 


When  the  dipoles  are  inclined  to  the  array  axis  at  the  special  angle 

sin2  0Q  - , (36) 

the  mutual  impedance  (34a)  simplifies  to  [12,14*] 

Zv«hW(kD  tel).  (37) 

The  active  impedance  is  then  given  by  the  closed  form 


«(£)  “ Jd~  W [n\}  [coskD  “ 1 

* — Indices  kD  -cosfclj,  l£l<». 


1 £ 1 < kD, 

(38a) 

9 . 

(38b) 

For  dipoles  spaced  more  than  1/2  wavelength  apart  (kD  > tr)  the  active  impedance  has  a 
nonvanishing  real  part  for  all  values  of  £.  All  surface-wave  roots  of  (17)  can  therefore  be 
found**  using  (38b): 


cos  * cos  kD  exp  (k£)X0) . (39) 

•Footnote  on  page  2. 

••SenguoU  ! 8 } derived  u aimilar  form.  However,  his  result  wu  intended  u an  approximation  for  root*  of  the 
conventional  Ys*i  structure,  with  Oq  * w/3. 
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In  the  preceding  section  (Eq.  (25))  the  general  solution  for  the  currents  I was  formu- 
lated as  a contour  integral  in  the  complex  f = plane:  n 


1 
2i Tj 


f. 


E Mn  + 1) 

j} 


}Xa  + 9(f) 


jdf 


Here 


(40a) 


in  wh'ch 


9(f)  =^Ln(l-of)  + l-^Ln(l-|) 


1 - 


kD 


Ln- 


o 


e-JkD  § 


(40b) 

(40c) 


(40d) 


with  9(f)  representing  the  continuation  of  (38)  into  the  complex  f plane,  the  contour  C 
bmng  whe  unit  circle  If  I*  l indented  as  shewn  in  Fig.  4a,  and  the  notation  Ln  f ' denoting 
the  principal  branch  of  tne  natural  ’ogarithm  function: 


Ln  f'  = ln  9'  + -ir  < < xr . 

The  active  impedance  function  9(f),  and  consequently  the  integrand  in  (40a),  has  four 
branch  points,  at 


f " 00,7,~ ,0, 


(42) 


The  branch  cuts,  corresponding  to  (41),  have  been  chosen  to  run  between  00  ami  a and 
between  l/o  and  0 The  integrand  also  his  poles  at  the  zeros  of  the  denominator 

j'(P)  n e , given  by  (39),  and  at  a po'e  ol  the  numerator  which  occurs  at  infinity  when 
n <-l  The  contour  C is  indented  about  the  branen  points  only  to  avoid  ambiguities  in  the 
drawing,  since  the  singularities  at  thos^  joints  are  integraLle.  On  the  other  hand,  tlu  hdenta- 
tions  at  the  sut  face- wave  poles  ff<p,are  essential  and  are  dictated  by  uniqueness  of  t ic  solu- 
tion for  the  currents.  A detailed  discuss  on  of  these  points  is  contained  in  trio  Appendix. 


Preparatory  to  evaluating  the  .ntegral  in  (40a),  the  contour  C may  be  deformed  about 
the  singularities  of  tire  integrand  as  shown  in  Fig.  4b  *oi  n > 0 and  as  shown  in  Fig.  4c  for 
n < 0.  In  either  case  the  integral  naturally  senates  into  the  two  components  given  by  (21a). 
A surface-wave  contribution  ltn  arises  front  the  pole,  and  a correction  contribution  / arises 
from  integration  around  the  branch  cut.  For  n > 0 the  circula*  arc  of  large  radius  in  S?ig.  4b 
docs  not  contribute. 
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The  general  form  of  the  contributions  due  to  the  (simple)  poles  has  already  been  stated 
as  (26a)  or  (26b).  Substitution  in  (26a)  of  the  active  impedance  function  (40c)  or  substitu- 
tion in  (26b)  of  the  active  impedance  (38b)  for  the  dipoles  inclined  at  the  special  angle  (36) 
yields 


Jsn  = -sgn(p)Eg{kD)e 


cos  kD  - cos  £ 


(P) 


sin  £ 


(P) 


(43) 


The  function 


sgn(p)  ■*  +1,  P>0, 

» -1,  P<0, 

accounts  for  the  sense  in  which  the  contour  circling  the  pole  is  traced.  In  view  of  the  relation 
(19)  which  is  verified  by  (39) 


‘sn 


s(-n) 


(44) 


as  expected  from  symmetry. 

The  integral  around  the  branch  cuts  must  be  evaluated  numerically.  For  n > 0 an  ap- 
propriate form  is  readily  obtained  by  use  of  (40b)  and  the  relations 

Ln(l-of)  ■ In  (r?  - 1)  + jit,  «>>t)>1,  (46a) 

for  the  inward  portion  of  the  path  of  integration  and 

Ln(l-of)51  In  (tj  - 1) (46b) 


for  the  outward  portion.  The  result  is 


For  n > 0 it  is  easy  to  estimate  the  variation  of  / . The  denominators  in  (46)  are  }X  + g{£) 
evaluated  near  the  branch  cut.  Each  denominator  has  some  finite  minimum  absolute  ^alue 
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independent  of  n.  It  does  not  take  on  the  value  zero;  that  value  occurs  only  at  the  roots 
Integration  then  yields  * 


cn 


f i£ 

2n  n 


n > 0, 


(47) 


where  A is  independent  of  n.  Thus  I Icn  I decreases  at  least  as  fast  as  1/n. 

For  n < 0 the  integrand  diverges  as  f -►  °°.  It  is  convenient  to  deform  the  contour  C as 
shown  in  Fig.  4c.  An  appropriate  form  of  the  branch-cut  integral  is  readily  formulated  by  use 
of  (40b)  and  the  relations 


Ln  ^1  ~j)  = In  (1  -17)  + ji r,  1 >v  > 0, 


for  the  inward  portion  of  the  path  of  integration  and 


Ln*1~f)  = In (1  — r?)  — 0<t?<1, 

for  the  outward  portion  of  the  path  of  integration.  The  result  is 
E Q+jnkD  rl  f «-n-l 


/ = ■ 

cn 


6+)>ikD  r'  t) 

2,1  •'o  fx-A)  * J*„  -1  i:k 


,-n-l 


l1  * w)  * ,x°  - ‘wu  [({  - *)  l1  - ■">)]  _ 


>dt). 


(48a) 


(48b) 


(49) 


Integration  then  yields 


cn 


Hi  8 
2tr  n ’ 


n <0, 


(60) 


where  B is  independent  of  n.  Thus  I Icn  I again  is  shown  to  decrease  at  least  as  fast  as  1/n.  In 
view  of  symmetry  ojily  one  of  the  integrals  in  (46)  and  (49)  need  be  evaluated.  In  fact,  using 
the  transformation  tj  = 1/tj,  one  verifies  directly  that 


!cn  m 'c(-n)  • (6D 

The  currents  / and  Icn  were  computed  for  a Yagi  structure  consisting  of  an  array  of 
dipole  antennas  inclined  to  the  array  axis  at  the  special  angle  On  ■ 64.74°  and  spaced 
D = 0.200  m apart.  Each  dipole  has  length  b =*  0.400  m and  radius  a = 0.00636  m.  Absolute 
lengths  are  specified  to  emphasize  that  the  curve  shows  the  true  frequency  dependence  of  a 
Yagi  structure  of  fixed  geometry.  Of  course  this  structure  can  l>e  scaled  in  the  usual  way 
from  the  assumed  frequency  or  design  wavelength  of  1 m.  The  dipole  radius  was  selected  to 
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allow  comparison  of  the  results  of  the  next  section,  for  = 90°,  with  those  of  Mailloux. 
The  normalized  self -reactance  of  the  dipole  X has  been  taken  from  Jordan  [19,20] . The 
input  reactance  is  normalized  by  dividing  by  the  input  resistance.  This  normalized  value  is 
the  same  whether  referred  to  the  loop  or  the  base  (center)  of  the  antenna.  The  formulas  for 
the  antenna  impedance  are  summarized  in  the  next  section  by  (53)  and  (54).  Significant 
current  ratios  are  plotted  in  Fig.  5 for  a range  of  free-space  wavelengths  around  1 m 
(feD  = 1.257). 

At  1 m or  less  the  self-reactance  of  the  dipole  is  high,  allowing  only  small  parasitic 
currents  to  be  induced  on  the  dipole  structure.  The  surface-wave  currents,  such  as  do  exist, 
propagate  at  nearly  free-space  velocity.  At  higher  frequencies  ( kD  > 1.257)  the  dipoles 
approach  resonant  length.  The  lower  self-reactance  allows  high  currents  to  be  induced  on  the 
dipoles.  The  surface  wave  is  strongly  coupled  to  the  slow-wave  structure  and  propagates  at  a 


Fig.  5 — Space-wave  (correction)  component  of  the 
current  relative  to  the  turfite-wave  component  at 
the  feed  point  and  the  relative  surface-wave  veloc- 
ity, computed  for  an  infinite  array  of  dipoles  inclined 
at  0O  ■ arctan  J~T  * 54.74°  with  dipole  radiua 
a • 0.00635  m,  dipole  length  b « 0.400  m,  and  di- 
pole spacing  D * 0.200  m 
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velocity  much  below  free-space  velocity.*  The  unattenuated  surface-wave  component  of  the 
current  becomes  large  relative  to  the  space-wave  component,  which  in  addition  decays  away 
from  the  excited  element.  The  rate  of  decay  of  the  space-wave  component  of  current  is 
shown  in  Fig.  6.  It  is  seen  to  be  slightly  more  rapid  than  the  lower  bound  1/n  developed  in 
(47)  and  (50). 


SURFACE-WAVE  CURRENTS  ON  GENERALIZED 
INFINITE  YAGI  STRUCTURES 

In  this  section  the  dispersion  characteristics  of  the  surface-wave  currents  on  arrays  of 
arbitrarily  oriented  dipoles  will  be  computed.  As  in  the  case  of  the  specially  inclined  dipoles 


Fig.  8 — Decay  along  the  Vagi  structure  of  the 
space-wave  component  of  the  current  relative 
to  the  surface-wavo  component,  computed  for 
an  infinite  array  of  inclined  dipoles  with  the 
name  0n , u,  b,  and  I)  as  in  Fig,  6 and  for  a free- 
space  wavelength  of  1 m (iti)>  1.257) 


♦The  characteristics  of  the  surface  wnu  are  also  found  by  interpolating  for  0q  ■ 54.74°  on  the  wavenumber 
diagram  to  be  presented  as  Fig.  7.  The  slope  of  a line  from  (lie  origin  to  a point  on  a dispersion  curve  is  the 
relative  phase  velocity  of  the  corresponding  surface  wave  v^/c,  given  by  ( 20). 
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of  the  preceding  section,  the  generally  oriented  dipoles  are  easily  covered  by  the  network 
formulation  given  in  the  second  section.  The  currents  excited  by  any  one  driven  dipole  are 
given  by  (15),  and  the  dispersion  relation  for  the  surface-wave  components  is  given  by  (17). 
Nevertheless  analytical  difficulties  arise  due  to  the  increased  complexity  of  the  active  im- 
pedance Q(|).  Although  these  difficulties  stand  in  the  way  of  a complete  solution  along  the 
lines  obtained  in  the  preceding  section,  the  dispersion  relation  is  readily  solved  numerically 
for  the  surface-wave  parameters. 

When  the  more  general  expression  for  the  mutual  impedance  (34a)  is  substituted  into 
the  formula  for  the  active  impedance  (14),  the  resulting  series  may  be  summed  by  an  exten- 
sion of  the  methods  employed  in  the  Appendix.  This  sum,  found  in  Ref.  12,  is 


ReQ®  * (^)(om280- 2^a<o)+ sw*taVtl<u,<t-  <52a) 

= 0,  kD<  1 1 1 < jr,  (52b) 


and 


Ira  Q($)  » u ({)  Jcos2  0Q  - j sin2  6^j  + w Q)  sin2  6Q  , 


(52c) 


where 


«(*>■ 


I fkD-t  rkD+i 

-^-—  2.404  + J ( kD  -£  - x)  In  2 Ein^dx  + J (kD  + £ -x)  In  2 sin^-dx 

. /.MM  rfcD-»  ] 

J ln2sin—dx  + J ln2sin~dx  , l£l<kD<ir, 


(wy 


(52d) 


(kD)* 


2.404  + 


'o 

/t  + kD  fl-kV 

ln2ain-|dx  - I ln2sin|-dx  , «Kl(l<», 
•'o 


/t-M> 

«- 


kD  - x)  In  2 sin 
\kD 


/(*kD 


+ kD  - x)  in  2 sin  ^dx 


(hoy 


(Me) 


'0  o 

and  (with  a factor  3/2  added  that  is  missing  in  Ref.  12) 

u>(£)  « -gj—  In  2 (cos £ - coikD),  \£\<kD<ir, 


2 kD 


In  2 (cos kD  - cos  $),  &D<l£l<t 


<ff. 

(52f) 

Kir. 

(52g) 
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The  reactance  of  the  dipole  k may  be  obtained  by  any  independent  means.  Herein 
results  quoted  from  Jordan  and  Balmain  [19,20]  are  used.  Within  the  approximations  em- 
ployed for  mutual  coupling,  the  physical  dimensions  of  the  dipole  (length  and  radius)  enter 
only  through  the  value  of  the  input  impedance  ka  + jkQ : 


and 


= 15[(2  + 2 cos  h)  S1(S)  - (cos  b)  ■S1(26) 

- 2(sin  S)  Si( S)  + (sin  b)  Si(2b)]  (53) 

jXa  = - 15  j(sin  6)  [ —7  + ln^  + 2Ci(S)  - Ci(2b)] 

- (cos  h [2  Si(b)  - Si(2h)}  - 2Si(S)l,  (54) 


where  b is  the  dipole  length  kb  in  electrical  radians  (b  being  the  dipole  length  in  meters),  3 is 
the  dipole  radius  ka  in  electrical  radians  (a  being  the  dipole  radius  in  meters),  a/b  < 1,  y is 
Euler’s  constant  (0.5772),  and  S^(x),  Ci(x),  and  Si(x)  are  defined  as 


(55a) 

(55b) 


and 


(55c) 


The  normalized  values  of  antenna  impedance 

*0  + #0  “ 1 + ^ A (66) 

do  not  depend  on  the  current  point  on  the  antenna  chosen  for  reference. 


The  characteristics  of  the  surface  wave  obtained  from  the  solution  of  the  dispersion  re* 
lation  (19)  may  bo  displayed  in  various  ways.  Perhaps  the  most  universal  way  is  a wavenumber 
diagram:  kD  vs  $t.  For  convenience  of  scale  in  Fig.  7 to  follow,  these  parameters  will  be 
divided  by  2v  to  obtain  D/K  and  D/Xg  (fractional  wavelengths). 
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The  general  features  of  the  wavenumber  diagram  are  well  known  [4] , and  special  aspects 
for  the  case  of  the  normal  Yagi  structure  (0Q  = 90°)  are  discussed  by  Mailloux  [21] . In  par- 
ticular Mailloux  discusses  the  possibility  of  wave  solutions  for  antennas  of  length  X/2  < b < X 
(inductive  antennas).  He  concludes  that  the  antennas  must  be  capacitive  to  support  a surface 
wave.  The  present  computations,  shown  in  Fig.  7,  agree  with  those  of  Mailloux  and  confirm 
this  conclusion  for  the  case  0^  = 90°.  The  region  of  capacitive  antenna  length,  approximately 
b A < 0.5,  is  restricted  to  a region  of  each  diagram  that  is  below 


if-  < 


(57) 


For  example,  in  Fig.  7a  this  capacitive  region  is  below  Dfk  = 0.5  D/b  - 0.5  (0.200)/0.300 
= 0.33. 


However,  when  0^  -■£  0,  it  is  evident  just  from  the  form  of  the  imaginary  part  of  Q(£), 
given  by  (52c),  that  inductive  dipoles  will  generally  lead  to  solutions  for  some  0..  less  than 
90°;  unless  u(£)  = 0,  Im  Q(£)  will  change  sign  as  -*■  0.  Figure  7 shows  wave  solutions  when 
bfk  > 0.5  for  0q  ~ 30°  and  0°.  In  the  computations  (Fig.  7)  the  mutual  coupling  between 
antennas  has  (apart  from  a scale  factor)  been  approximated  by  the  mutual  impedance 
between  eanonical-ininimum-seattering  antennas  having  the  radiation  pattern  (in  the  open- 
circuited-array  environment)  of  short  dipoles  (Kqs.  (34)  and  (52)).  Therefore  results 
obtained  for  inductive  dipoles  apply,  strictly  speaking,  only  to  short  inductively  loaded 
dipoles.  Because  of  the  slow  change  of  the  radiation  pattern  with  dipole  length,  however, 
the  theory  of  mutual  coupling  between  minimum-scattering  antennas  [17]  indicates  that 
these  results  are  approximately  applicable  to  arrays  of  (unloaded)  dipoles  somewhat  longer 
than  1/2  wavelength.  In  Fig.  7 the  input  impedance  given  by  (53)  and  (54)  is  specified  in 
terms  of  a fixed  dipole  length  6,  fixed  dipole  radius  a,  and  the  free-spaee  wavelength  X.  Of 
course  the  form  of  (54)  plays  a large  role  hr  determining  the  frequency  dependence  shown. 

A part  of  the  active  impedance  He  Q(£)  in  the  range  0 < 1 $ I < X7>,  although  it  does 
not  enter  into  the  resonance  calculations  (Re  Q(£)  s 0.  t:D  < I £ I < tt),  is  readily  compared 
for  short  and  1/2-wavelength  dipoles.  This  portion  of  the  active-impedance  summation  can 
be  rigorously  expressed  in  terms  of  the  radiation  patterns  [12] 


cos 


where  RQ  is  the  input  resistance  to  the  antenna  element  in  the  open-circuited  array  environ-, 
ment.  For  a CMS  antenna  /i0  =*  li , which  has  been  normalized  to  1.  In  (58)  Pa{Q,$)  is  the 
radiation  pattern  of  an  element  in  me  open-circuited-array  environment,  normalized  to  unit 
radiated  power.  For  the  readily  computed  special  case  of  coliinear  dipoles  (0Q  ■ 0) 

[‘  -(41  l59) 


for  short  dipoles  (52a)  and 
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Re  Q(£) 


(60) 


for  half-wave  dipoles  [ 12,17,19] . (In  (60)  an  error  in  the  coefficient  of  the  corresponding 
expression  (37)  of  Ref.  12  and  (78)  of  Ref.  17  has  been  corrected.)  Although  the  functional 
form  is  different,  numerically  the  patterns  and  consequent  variations  of  Re  Q(£)  are  similar. 
In  the  range  0 < l £ I < 0.6  kb  the  two  differ  by  less  than  10%. 


Figure  8 illustrates  several  features  of  significance  in  the  computations.  The  first  feature 
is  ihe  dependence  of  the  active  reactance  Im  Q(£)  on  £.  The  second  is  the  solution  of  the 
resonance  relation  £s  at  the  intersection  of  the  negative-antenna-reactance  {-Xa)  line  with 
the  active-reactance  line.  The  third  is  the  extreme  sensitivity  of  the  solution  (£s)  to  small 
changes  in  antenna  reactance  for  larger  values  of  £ (the  slowest  waves).  Because  of  this 
sensitivity  the  solution  becomes  unreliable;  hence  the  dispersion  curves  in  Fig.  7 have  been 
broken  off. 


Independent  results  with  which  the  present  theory  may  be  compared  are  available  for 
only  the  conventional  Yagi  structure,  with  dQ  = 90°,  and  consequently  only  for  antenna 
reactances  corresponding  to  inductive  dipole  lengths  b < A/2.  A comparison  with  the  theory 
of  Mailloux  [7]  is  shown  in  Fig.  9.  The  relatively  simple  expression  for  antenna  reactance 
used  (Eq.  (54))  is  at  about  the  limit  of  its  range  of  validity  in  terms  of  dipole  radius  [19] . 

A possible  application  of  the  new  information  on  Yagi  structures  contained  in  Fig.  7 is 
in  design  of  a mechanically  compensated  broadband  Yagi  antenna.  The  well  known  condition 
for  optimum  directivity  of  a Yagi  surface-wave  antenna  (Hansen-Woodyard)  is  cited  by 
Ehrcnspeck  and  Poehler  (3) . In  the  notation  of  Fig.  7,  the  required  relation  between  wave- 
lengths and  the  length  L of  the  Yagi  director  array  is 


0.468 


2 *D 

r 


D 


(61) 


In  Pig.  10  this  relation  is  shown  os  the  curve  labeled  “optimum"  superimposed  on  a (schem- 
atic) wavenumber  diagram  for  a generalised  Yagi  structure.  As  frequency  (DA)  is  increased, 
a conventional  director  structure  ((?Q  « 90°)  of  fixed  dipoles  proportioned  to  satisfy  relation 
(61)  (point  A ) rapidly  departs  from  this  condition  to  point  B.  However,  if  each  dipole  is 
rotated  appropriately  to  an  inclination  ft,  < 90° , the  optimum  relation  is  restored  at  the 
higher  frequency  (point  C).  Thus  a director  array  of  dipoles,  mechanically  ganged  so  that 
frequency  and  inclination  track,  produces  optimum  performance  of  the  director  at  each 
frequency.  However,  the  concomitant  rotation  of  the  element  pattern  for  each  dipole  does 
involve  some  loss  of  directivity. 
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Fig.  8.  — Calculated  effects  of  the  surface-wave  phase-shift- 
parameter  root  £s  for  a conventional  Yagi  structure  (0o  ■ 90  ) 
with  a = 0.00635  m,  b = 0.440  m,  D • 0.200  m.  and  O/X  = 0.196 
(feD-  1.23) 
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Fig,  9 — Comparison  of 
the  relative  phase  veloc- 
ity of  a surface  v,  -re 
from  the  theory  of  Mali- 
loux  [7]  (solid  curve) 
and  from  the  present 
analysis  (points) 


Fig.  10  — Rotation  of  the  dipole  inclination  Oci 
required  for  frequency  (7)/X)  compensation  of 
a Yagi  director  array.  When  $0  is  90°  and  per- 
formance ia  optimum  for  a first  frequency 
(point  A),  the  performance  la  not  optimum  for 
a higher  frequency  (point  B).  However,  decreas- 
ing the  dipole  inclination  to  an  angle  0o  < 90° 
would  again  give  optimum  performance  (point 
C)  at  the  higher  frequency. 
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APPENDIX 

The  Functional  Forms  of  Q(£)  and  q(f) 

Symmetry  and  Lorentz  reciprocity  result  in  relation  (4)  among  the  impedance 
coefficients: 


*,  m Z-v  ’ 


As  a result  the  active  impedance  as  given  by  (14)  becomes 


(Al) 


W) 

Consequently,  if  the  equation 


JX0  * <?(*)  - 0 

has  a root  , then  it  has  another  root  such  that 

t<-P>  • 

's  S 

This  is  a general  feature  of  periodic  reciprocal  structures  14,22] . 


(A2) 

(A3) 


(A4) 


25 


KAHN 


If  the  complex  variable  £ = ije^  is  introduced,  the  complex  active  impedance  function 


«(f)  - £ zvrv 

v-— 

can  be  put  into  the  form 

4(f)  ■ zo*  tl  ZX  * r'j. 

v-1 

Consequently,  if  the  equation 


jXa  + q(f)  = 0 

has  a root  then  it  has  another  root  such  that 
s * 


(A5a) 


(A5b) 


(A3) 


►(-P)  - 1 
S - f(P)  * 


(A7) 


The  active  impedance  (A2)  can  be  summed  in  closed  form  for  an  array  of  dipoles 
inclined  at  the  special  angle  sin2  ■ 2/3,  at  which  angle  Z reduces  to  (37),  which  by  use 
of  (35a)  is 


e-]\v\kD 

Z'>  " ; \»\kD  ‘ 


(A8a) 


In  the  case  v » 0 

ZQ  « 1. 


(A8b) 


The  remainder  of  this  appendix  deals  only  with  this  special  anglo.  Substituting  (A8a)  and 
(A8b)  in  (A2)»  one  has 


Qtt) 


.-MkD  -|) 


+ 1 + 


4r% 


,-MkDH) 


If  k is  assumed  to  have  a small,  negative  imaginary  part,  so  that 

* " *r  + kt  < 

then  the  scries 


(A9) 


(A10) 


i>“0  »>-0 


(AH) 
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is  absolutely  convergent.  Integrating  term  by  term  (limits  0 and  of)  yields  the  first  summa- 
tion in  (A9).  On  the  other  hand  (All)  is  a geometric  progression  with  the  sum 

■nhr-  <Al2> 

Consequently  the  first  sum  in  (A9)  is  also  given  by  the  integral  of  (A12): 


Yzhc  = — ln  (!  — (A13) 

in  which  In  f ' denotes  the  natural  logarithm  of  £'.  The  expression  (A13)  for  the  sum  is  a 
multivalued  function  of  of,  since  In  f ' is  multivalued.  But  the  sum  (All)  is  single  valued;  it 
converges  to  the  principal  value  of  the  logarithmic  function  Ln  f ' defined  by 

-rr  < Im  (Ln  f ')  < +7r.  (A15) 

The  second  summation  in  (A9)  may  be  performed  similarly.  The  complete  result  is 

= ~kD  Ln(1  ~af)  + 1 ~~kDln(1  ~f)-  (A16) 


In  the  limit  k(  -*■  0 


a(f)  » » Q(£).  (A17) 

In  the  main  text  the  finite  Fourier  transform  was  employed  to  obtain  a formal  solution 
of  the  difference  equation  (5)  for  the  (Green’s  function)  currents  /„ . The  result  is 


i r 

2*4  iXa  ♦ Q(f) 


(A18) 


where  the  value  of  fi  in  (15)  has  been  set  equal  to  zero  and  this  superscript  suppressed.  This 
solution  was  obttiined  without  explicit  reference  to  boundary  conditions  (in  this  case  condi- 
tions at  I n I -+  <*> .)  Without  such  boundary  conditions  or  their  equivalent,  any  solution  of 
the  homogeneous  difference  equation  may  be  added  to  a given  solution  for  the  currents, 
producing  another  solution.  This  lack  of  uniqueness  manifests  itself  through  the  presence 
of  singularities  on  the  path  of  integration.  One  might  intuitively  associate  a zero  of  the 
denominator  with  finite  /„  even  for  E.  -*•  0,  that  is,  with  solutions  of  the  homogeneous 
equation. 

A unique  solution  is  obtained  when  nonzero  solutions  of  the  homogeneous  equations 
are  excluded  by  introducing  dissipation.  Some  dissipation,  no  matter  how  little,  is  physically 
unavoidable;  therefore  this  process  selects  the  physically  correct  unique  solution  in  the 
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limit  as  the  dissipation  approaches  zero.  In  terms  of  the  singularities  of  (A18),  dissipation 
removes  these  singularities  from  the  real  path  of  integration  -it  < £ < it  to  complex  values. 
The  movement  of  the  singularities  will  now  be  studied  in  the  complex  f plane. 

The  singularities  in  the  integrand  of  ( 40a)  are  the  branch  points  of  g(f):  roots  of 
jXa  + q($)  = 0 and  singularities  of  f“(n+r).  The  branch  points  of  the  function  qr(f)  arise 
from  branch  points  of  Ln  (f')  in  (A16).  These  are  listed  in  Table  Al.  As  is  well  known,  dis- 
sipation produces  a negative  imaginary  component  in  k,  so  that  (A10)  applies.  The  effect  is 
to  move  the  branch  points  1 la  and  a off  the  unit  circle: 


and 


il/ol  = \ejkD  1 = e“*i'-*  > 1 

(A19a) 

1 a 1 = \e~ikD  1 = < 1. 

(A19b) 

Table  Al  — Branch  Points  of  the  Active 
Impedance  Function  q(f ) 


Branch  Points  of  Ln  f ' 

Branch  Points  of  q (f) 

f'  = 1 - of  = 0 

f*  **  1 — of  -*■  oo 

f'  = 1 - o/f  ■ 0 

f'  = 1 - o/f  -*  “> 

f - l/o 

branch  cut 

f 00 

f = a | 

> branch  cut 
f =>  0 l 

Properties  of  the  roots  of  (A3)  and  (A4)  have  already  been  discussed.  Displacement  of 
the  roots  is  readily  calculated  when  dissipation  is  introduced  in  the  form  of  a (small)  real 
component  in  the  antenna  self-impedance.  Then  (A4)  becomes 

Ra  + iXa  * 9(f)  - 0.  (A20) 


Expanding  g(f)  in  a Taylor  series  about  an  unperturbed  root  f results  in 


9(f) 


(A21) 
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By  construction,  to  first  order  in  f - f ^ , one  has 

s 


or 


R‘  * f M [f 

y-(P) 

f - if*  - -«.( 

_1 

dtl 

= 0, 


(p) 


(A22a) 


(A22b) 


The  right-hand  side  of  (A22b)  has  been  evaluated  in  another  connection;  from  (26a)  and 
(43),  with  n = -1,  the  result  is 


In  (A23),  forp  = +1 

and  for  p = -1 


r - 


if’  - *. 


kDt 


(P) 


cos  kD  - cosf 


(P) 


sin  f 


(P) 


feD<C(1)  < 7T 
$ 

-ir  < t(-1)  < kD. 


(A23) 


(A24a) 

(A24b) 


The  trigonometric  factor  is  positive  in  the  case  of  (24a)  and  negative  in  the  case  of  (24b). 
Consequently  the  root  H1*  moves  radially  outward  (along  H1))  and  the  root  moves 
radially  inward  (along  -K”*)).  The  location  of  the  singularities  in  the  presence  of  dissipa- 
tion is  illustrated  in  Fig.  Al. 


Fig.  Al  — Singularities  of  the  integrand  in  (40a)  in  the 
complex  f plane  in  the  pretence  of  dluipation  when 
tlnJfl0  - 2/3 
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